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Abstract 


YTwo  theoretical  descriptions  have  been  developed  for  the  phase  boundary 
dynamics  during  crystal  1 1zatlon  of  amorphous  films  by  scanning  with  the  slit 
image  of  a  cw  laser  or  of  any  cw  energy  beam.  The  first  reduces  the  problem 
to  the  solution  of  a  one-dimensional  Integral  equation,  which  allows  a  choice 
of  Initial  conditions.  Depending  on  the  background  temperature,  numerical 
solutions  yield  either  periodic  or  runaway  motion  of  the  amorphous-crystalline 
(a-c)  boundary,  as  observed  In  experiments  on  scanned  laser  crystallization  of 
thin  films  of  a-Ge  on  fused-slllca  substrates.  The  calculations  give  a 
seml-quantltatlve  fit  to  the  experimental  results  for  the  spatial  periodicity 
observed  In  the  crystallized  films  as  a  function  of  background  temperature. 
Profiles  of  film  temperature  as  a  function  of  distance  from  the  laser  Image  at 
successive  times  have  been  computed  for  both  the  periodic  and  runaway  cases. 
The  model  qualitatively  explains  many  of  the  effects  observed  during  scanned 
cw  laser  crystallization.  Including  periodic  fluctuations  In  light  emission. 
The  second  theoretical  description  Is  a  more  exact  two-dimensional  treatment, 
applicable  only  to  cases  of  steady-state  motion  of  the  a-c  boundary,  which 
rigorously  handles  heat  flow  Into  the  substrate.  This  treatment  has  been  used 
to  calculate  the  boundary  velocity  during  steady-state  runaway .^The 
dependence  of  this  velocity  on  background  temperature  and  on  film  and 
substrate  thermal  properties  and  thickness  has  been  determined.  At  the 
minimum  background  temperature  required  for  runaway  the  calculated  value  of 
the  steady-state  velocity  is  ~  140  cm/sec  for  the  case  of  a  Ge  film  0.3  un 
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thick  on  a  fused-sillca  substrate  1  mm  thick.  Experimental  values  lie  In  the 
range  100-300  cm/sec.  A  class  of  laser-controlled  steady-state  solutions  has 
been  obtained  for  which  the  boundary  velocity  Is  equal  to  the  laser  scanning 
velocity  but  lower  than  the  boundary  velocity  for  uncontrolled  runaway.  The 
existence  of  these  solutions  suggests  the  possibility  of  preparing 
single-crystal  semiconductor  sheets  by  scanning  amorphous  films  In  a  manner 
that  achieves  uniform,  laser-controlled  motion  of  the  a-c  boundary  rather  than 
either  periodic  or  runaway  motion. 
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ILLUSTRATIONS 


Schematic  representation  of  laser  crystallization  experiment. 

Normalized  position  of  a-c  phase  boundary,  S(t),  as  a  function  of 
normalized  time,  t,  given  by  solutions  of  Eq.  (3)  for  three 
different  values  of  n.  The  Inset  shows  the  result  of  an 
experiment  In  which  change  of  transmission  is  proportional  to 
S(  t). 

Optical  transmission  micrograph  of  a  laser-crystal llzed  Ge  film 
showing  periodic  structural  features.  Inset  Is  expanded  view 
illustrating  four  different  microstructure  regions. 

Left:  Brlght-field  transmission-electron  micrograph  of  a 
laser-crystallized  Ge  film.  Illustrating  four  different 
microstructure  regions.  Right:  Transmission-electron 
diffraction  patterns  for  fine-grained  and  large-grained  regions 
(upper  and  lower,  respectively). 

Vat  vs  n  from  solutions  of  Eq.  (3). 

Theoretical  curve  and  experimental  points  for  AY/AY0  vs 
temperature  ratio  Tjj/Tr.  Curve  Is  calculated  using  Fig.  5  with 
no  »  0.22  and  Tc  *  680°C. 

Normalized  temperature  distribution  T'(S',  t)  as  a  function  of  S' 
for  several  values  of  t.  Parameters  were  chosen  to  Illustrate 
the  case  of  periodic  a-c  boundary  motion. 

Normalized  temperature  distribution  T'(S',  t)  as  a  function  of  S' 
for  several  values  of  t.  Illustrating  periodic  temperature 
flare-up. 

Normalized  temperature  distribution  T'(S',  t)  as  a  function  of  S' 
for  several  values  of  t.  The  parameter  values  o*  2.0,  n*  0.8, 

V  ■  0.3  and  r  »  30  have  been  chosen  to  Illustrate  the  case  of 
runaway. 

F(S ',  v)  as  a  function  of  S'  for  several  values  of  V. 
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11  Plots  of  log  [F(0.  V)]  vs  log  V  for  6  *  240,  120,  and  60,  and 
plots  of  log  [n1'2  n*V]  vs  loq  V  for  n'  *  9.0  and  n'  *  2.2. 

12  Blown  up  view  of  periodic  features  In  a  Ge  film  crystallized  by 
the  slit  Image  of  a  laser.  Picture  takei  near  the  center  axis  of 
the  slit  Image  shows  the  origins  of  chevron  pattern. 

B1  Diagram  for  the  calculation  of  the  exact  stationary  solution  of 
moving  a-c  boundary  problem. 


I.  Introduction 

The  study  of  the  transformation  of  semiconductor  films  from  the 
amorphous  to  the  crystalline  state  has  become  a  matter  of  great  Interest  In 
the  last  few  years.  The  time  dependence  of  the  transformation  has  been 
investigated  as  a  function  of  film  temperature  for  Ge1  and  SI ,2  and  the  latent 
heat  of  transformation  and  transformation  temperature  has  been  measured  for 
Ge,3*4  SI,4  and  Ge-Si  alloys.4  The  so-called  "explosive"  transformation  of 
semiconductor  films  has  been  reported  by  a  number  of  workers, 5  and  the 
velocity  of  the  transformation  has  been  measured. 6 

We  have  recently  reported  the  observation  of  a  number  of  unusual 
phenomena  during  a  study  of  the  laser  crystallization  of  amorphous  Ge  films.7 
These  phenomena  Included  the  formation  of  periodic  structural  features, 
pulsations  of  film  temperature  during  laser  scanning,  and  runaway 
crystallization  of  the  entire  film  following  momentary  contact  with  the  laser 
Image.  To  provide  a  qualitative  description  of  the  laser-crystallization 
process,  we  presented  a  one-dimensional  integral  equation  description  of  the 
amorphous-crystalline  (a-c)  phase  boundary  motion,8  which  takes  Into  account 
the  latent  heat  emitted  during  the  a-c  transformation.  Numerical  solution  of 
the  Integral  equation  provides  considerable  Insight  into  the  dynamics  of  the 
a-c  boundary.  The  results  were  used  to  obtain  a  semlquantltatlve  fit  of  the 
temperature  dependence  of  the  structural  periodicity  observed  In 
crystallization  of  amorphous  Ge  films  on  fused  silica  substrates. 
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In  this  report  we  present  a  somewhat  Improved  version  of  the 
one-dimensional  Integral  equation  description  of  laser  crystallization  which 
takes  account  of  the  loss  of  heat  to  the  film  substrate  by  Introducing  Into 
the  equation  an  ad-hoc  exponential  time  dependent  damping  factor.  The  results 
obtained  for  the  motion  of  the  a-c  boundary  are  qualitatively  very  similar  to 
those  obtained  earlier.®  The  solution  of  the  Integral  equation  Is  used  to 
obtain  series  of  plots  of  temperature  as  a  function  of  lateral  position  in  the 
film  at  fixed  times  after  the  onset  of  the  crystallization  process.  These 
plots  provide  an  understanding  of  the  temperature  pulsations  during  the 
formation  of  periodic  structures  In  film  morphology,  as  well  as  a  vivid 
picture  of  the  onset  of  a-c  phase  boundary  runaway.  The  general  features  of 
the  model  are  strongly  supported  by  our  recent  experimental  study  of  a-c 
boundary  dynamics  during  laser  crystallization.® 

The  one-dimensional  Integral  equation  solutions  provide  an  approximate 
description  of  transient  as  well  as  steady-state  behavior  of  the  a-c  phase 
boundary  motion.  On  the  other  hand,  we  have  obtained  a  class  of  less 
approximate  two-dimensional  solutions  which  are  stationary  in  a  frame  of 
reference  moving  with  the  a-c  phase  bounary.  It  Is  far  more  complicated  to 
try  to  obtain  transient  solutions  of  the  two-dimensional  problem.  The 
stationary  solutions  are  exact  If  the  temperature  dependence  of  film  and 
substrate  properties  as  well  as  differences  In  amorphous  and  crystalline 
properties  are  neglected.  The  solutions  allow  a  calculation  of  the  a-c 
runaway  boundary  velocity,  as  well  as  the  velocity  of  the  boundary  when  the 


2 


, 

.  *N.’  .  A 


laser  controls  its  motion.  The  runaway  velocities  calculated  are  In 
reasonable  agreement  with  observation. 

In  Section  II  of  this  paper  we  present  the  transient  and  periodic 
steady-state  solutions  of  the  one-dimensional  Integral  equation  description  of 
laser  crystallization.  The  model  results  are  compared  to  experimental  results 
for  laser  crystallization  of  a-Ge.  In  Section  III  we  obtain  and  compare  the 
steady  state  runaway  solutions  of  the  Integral  equation  approximation  and  the 
two  dimensional  model.  In  section  IV  we  consider  steady  state  laser-guided 
solutions  obtained  from  the  Integral  equation  approximation  and  the  exact 
model.  In  Section  V  the  results  obtained  and  the  possibility  of  achieving 
steady  state  laser-guided  a-c  boundary  motion  are  discussed.  The  numerical 
solution  of  the  Integral  equation  is  discussed  in  Appendix  A,  and  the 
stationary  solution  of  the  exact  two-dimensional  problem  is  derived  in 
Appendix  B. 

II.  Transient  and  Periodic  Solutions  of  the  Integral  Equation 

A.  Derivation  of  the  Integral  Equation 

When  the  temperature  of  an  amorphous  semiconductor  film  Is  raised, 
tranformatlon  to  the  stable  crystalline  form  takes  place  at  a  rate  that 
Increases  exponentially  with  temperature,* *2  so  that  over  a  narrow  temperature 
Interval  at  a  temperature  ~  Tc  the  ratio  of  the  time  taken  for  a  laser  scan 
changes  from  »1  to  «1.  We  therefore  argue  that  a  reasonable  description  of 
the  tranformatlon  Is  given  by  assuming  that  It  occurs  when  the  amorphous  film 
reaches  a  critical  temperature  Tc.  Support  for  this  point  of  view  comes  from 
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the  observation  of  a  rather  sharply  defined  transformation  temperature  in 
latent  heat  measurements. 4  it  has  been  suggested*0  that  the  transformation  may 
In  fact  correspond  to  the  change  from  the  amorphous  to  the  liquid  state, 
followed  by  a  transformation  to  the  crystalline  state.  In  this  case  the 
existence  of  a  well  defined  transition  temperature  Tc  is  expected.  The 
assumption  of  a  transformation  at  a  constant  temperature  Tc  should  be  a 
reasonable  approximation  In  either  case,  and  Is  justified  by  the  good 
description  of  experiment  that  results.  The  modification  of  this  assumption, 
and  the  detailed  nature  of  the  transformation  are  considered  further  in  the 
discussion  of  Sec.  V. 

We  assume  that  the  temperature  dependence  of  the  film  properties  and  the 
difference  In  these  properties  between  the  amorphous  and  crystalline  states 
can  be  neglected.  This  Is  an  acceptable  approximation  at  temperatures  In  the 
range  of  the  amorphous-crystalline  transformation  (~  500-700°C  for  Ge).  With 
these  simplifying  assumptions,  an  integral  equation  formalism*!  can  be  used  to 
describe  the  motion  of  the  phase  boundary.  The  geometry  assumed  for  the  laser 
crystallization  calculation  Is  shown  schematically  in  Fig.  I.  The 
semiconductor  film,  which  Is  deposited  on  a  thick  substrate,  is  of  Infinite 
extent  In  the  y  and  z  directions  and  so  thin  that  Its  temperature  Is  constant 
In  the  x  direction.  The  laser  slit  Image  Is  of  Infinite  length  In  the  z 
direction  and  moves  at  a  velocity  v  In  the  positive  y  direction.  At  time  t  = 

0,  the  phase  boundary  Is  located  at  y0,  with  the  crystalline  phase  to  the  left 
(y  <  yo)  and  the  untransformed  amorphous  phase  to  the  right  (y  >  y0).  The 
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Fig.  1.  Schematic  representation  of  laser  crystallization  experiment. 
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laser  image  carries  with  It  a  steady  state  temperature  profile  T0  (y-vt).  At 
t  =  0,  the  temperature  at  the  phaseJ>oundary  reaches  Tc  and  the  boundary 
begins  to  move  Irreversibly  toward  the  right,  with  heat  being  liberated  at  a 
rate  per  boundary  unit  cross  sectional  area  of  fLpV(t),  where  L  Is  the  latent 
heat  of  a-c  transformation  of  the  semiconductor,  p  Is  the  semiconductor 
density,  Y(t)  Is  the  position  of  the  boundary  at  time  t,  and  we  have  Included 
a  factor  f  which  is  less  than  1  and  accounts  In  an  approximate  way  for  the 
loss  of  a  fraction  (1-f)  of  the  latent  heat  to  the  substrate  as  it  Is 
liberated.  In  the  present  version,  we  will  also  Include  an  exponential 
damping  factor  y  to  represent  the  eventual  diffusion  Into  the  substrate  of  the 
fraction  f  of  the  heat  propagating  In  the  film.  We  shall  see  that  f  and  y, 
which  are  Introduced  as  ad  hoc  parameters  In  the  one-dimensional  Integral 
equation,  have  their  counterparts  In  the  two  dimensional  steady  state  model 
(Sec.  1 1 IB) ,  where  they  are  determined  In  the  course  of  the  calculation, 
leaving  no  adjustable  parameters.  The  effect  of  strain  and  other  mechanical 
forces  Is  not  explicitly  Included,  although  L  could  Include  a  contribution  due 
to  strain.  Then  the  temperature  T(y,t)  at  any  point  y  along  the  film  at  time 
t  Is  given  by  the  one-dimensional  Integral  relational 


T(y,t)  =  T0  (y-vt)  + 

t  t  r  i 

fQ  Y(t')  exp  {- y  [t-t']> 


(1) 


jexp  [4xx3"1/2  [t-tTl/2  dt' 


where  C  Is  the  specific  heat  of  the  film,  K  -  K/Cp  defines  Its  thermal 
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dlffuslvlty,  K  Is  Its  thermal  conductivity,  and  y  Is  a  phenomenological 
damping  factor. 

Equation  (1)  has  a  simple  physical  Interpretation.  It  states  that  the 
temperature  at  a  point  y  at  time  t  Is  a  superposition  of  the  contribution  due 
to  the  moving  laser,  T0  (y-vt),  and  the  sum  of  contributions  due  to  sources  of 
heat  fLpV(t')dt',  emitted  at  positions  Y(t')  at  earlier  times  t'.  The  source 


function  or  Green's  function 


[y  -  *(t'n 

4  it  -  t  ']  j 


-1/2  -1/2 
[4vk]  [t-t "] 


describes  the  one-dimensional  diffusion  of  heat  away  from  the  source,  and  the 
factor  exp  |-r[t-t']|  represents  the  decay  of  heat  out  of  the  film  and  into  the 
substrate  or  atmosphere. 

An  Integral  equation  for  Y(t)  can  be  obtained  by  using  the  condition 
that  the  temperaure  at  the  phase  boundary  Is  Tc,  or 

T  [Y(t),t]  -  Tc  .  (2) 


For  purposes  of  calculation.  It  Is  convenient  to  rewrite  Eq.  (2)  In  the  frame 
of  reference  moving  with  the  laser  Image.  We  Introduce  the  position  variable 
u(t)  *  y(t)  -vt,  where  u(t)  Is  measured  from  the  center  of  the  laser  Image  as 
origin.  The  temperature  T0(u)  Is  modeled  In  the  form  T0(u)  *  Tb  +  ATt  exp 
[-(u/a)2],  where  Tb  Is  a  uniform,  time- Independent  background  temperature  and 
the  temperature  contribution  due  to  the  laser  Is  described  by  a  Gaussian  of 
width  a  and  magnitude  AT Finally,  Introducing  a  unit  of  length  d  for 
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normalization  purposes  and  normalizing  quantities,  we  can  write  the  Integral 
equation  for  the  motion  of  the  phase  boundary  as 


1  *  “exp  J.  [SCt)]2| 


+  nf  C$(  t-)  +  V]/(t-t-)1/2}  x 


t'))  e 


exp  {-r(t-T'HLixP  {-[S (  t)-S( t')  +  V  (t-t')]2  /  (t-t')»t' 


where  a  *  AJV(TC-Tb),  (4a) 

h  ■  fn't  n'  *  L/Cirl/2(Tc-Tb) ,  (4b) 

T  =  4  <t/d2,  V  =  dv/4 ie,  S  s  U/d,  §  (t')  s  [dS( t')/<1t'],  r  «  (d2/4,c)y,  b  -  a/d, 
and  U  Is  the  position  of  the  phase  boundary  measured  from  the  center  of  the 
laser  Image.  For  the  discussion  of  this  section.  It  Is  convenient  to  choose 
d  =  a,  so  that  b  *  1. 

Equation  (3)  Is  not  In  Itself  sufficient  to  give  physically  acceptable 
solutions  for  the  motion  of  the  phase  boundary,  since  It  allows  negative  values 
of  [§(t')  +  V],  which  Imply  the  unphysical  motion  of  the  phase  boundary  back 
toward  the  laser  Image,  with  the  reconversion  of  crystalline  material  to  the 
amorphous  state,  accompanied  by  the  reabsorption  of  latent  heat.  To  constrain 
Eq.  (3)  to  physically  acceptable  solutions,  we  require  that  when  the  numerical 
solution  of  Eq.  (3)  yields  [§(t')  +  V]  <  0,  this  quantity  Is  to  be  set  equal 
to  zero,  with  the  phase  boundary  remaining  stationary. 
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Equation  (3)  has  been  solved  numerically  (see  appendix  A)  to  obtain  S  as 


a  function  of  t  for  representative  values  of  a,  n»  V,  and  r,  with  the  boundary 
conditions  that  §(t)  =  0  at  t  «  0  and  that  S(0)  Is  given  by  1  *  o  exp 
{-[S(0)]2}.  The  values  of  the  parameters  have  been  chosen  for  convenience  In 
numerical  solution  of  Eq.  (3)  rather  than  an  optimal  fit  to  experimental 
conditions.  Figure  2  shows  plots  of  S  vs  t  for  V  =  0.3,  r  *  30  and  three 
increasing  values  of  n,  with  a  Increasing  In  proportion  to  n  (which 
corresponds,  from  Eqs.  (4a)  and  (4b),  to  Increasing  Tj>  toward  Tc  while  holding 
AT*  fixed.)  For  each  value  of  n,  S  initially  Increases  rapidly  because  the 
latent  heat  liberated  by  the  phase  transformation  raises  the  temperature  ahead 
of  the  boundary,  accelerating  Its  forward  motion.  As  the  boundary  moves  away 
from  the  laser  Image,  the  contribution  of  the  laser  to  the  temperature  ahead  of 
the  boundary  decreases  rapidly.  For  n  *  0.3,  the  boundary  motion  soon 
decelerates,  and  the  boundary  comes  to  rest,  remaining  fixed  for  a  time 
Interval  during  which  Its  temperature  begins  to  drop  rapidly  below  Tc  and  S 
decreases  with  velocity  §  ”  -V.  With  the  approach  of  the  laser  Image  the 
boundary  temperature  gradually  Increases  to  Tc,  the  boundary  once  more  moves 
forward,  and  S  again  Increases.  This  cycle  Is  repeated  indefinitely, 
resulting  In  the  oscillations  In  S  seen  for  n  *  0.3  In  Fig.  2.  The  Inset  In 
Fig.  2  shows  the  result  of  an  experiment 9  which  measures  the  Infrared  light 
transmitted  through  a  thin  film  of  Ge  as  the  a-c  boundary  moves  In  a  cw  laser 
crystallization  experiment.  Since  the  transmission  of  crystalline  Ge  Is  much 
greater  than  that  of  amorphous  Ge,  the  changes  In  transmitted  light  signal 


9 


C83-4702  1 


Fig.  2.  Normalized  position  of  a-c  phase  boundary,  S(t),  as  a  function 
of  normalized  time,  x,  given  by  solutions  of  Eq.  (3)  for  three  different 
values  of  n.  The  Inset  shows  the  result  of  an  experiment  In  which  change 
of  transmission  Is  proportional  to  S(x). 
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with  time  are  proportional  to  the  a-c  boundary  position  with  time.  The 
similarity  between  the  oscillations  observed  In  the  experiment  and  those  shown 
In  the  curve  for  n  *  0.3  Is  evident. 

When  n  Is  Increased  by  Increasing  T^,  less  heat  Is  required  to  raise 
the  temperature  of  the  film  to  Tc  ahead  of  the  laser  Image,  and  the  phase 
boundary  moves  farther  beyond  the  laser  before  It  decelerates  and  comes  to 
rest.  This  trend  Is  Illustrated  by  the  calculated  curve  for  n  3  0.6  In  Fig.  2 
and  leads  to  motion  with  a  longer  period.  When  n  becomes  large  enough,  the 
heat  liberated  during  crystallization  is  sufficient  to  sustain  the 
transformation,  causing  the  boundary  to  "run  away"  from  the  laser  Image.  This 
situation,  which  Is  Illustrated  by  the  curve  for  n  ■  0.8  In  Fig.  2,  accounts 
for  the  observation  that  for  high  enough  background  temperature  the  entire 
film  Is  crystallized  following  momentary  contact  with  the  laser  Image. 

In  addition  to  explaining  runaway  crystallization,  the  proposed  model 
can  also  explain  our  other  qualitative  observations  on  laser  crystallization. 
The  model  does  not  directly  predict  observable  structural  changes  In  the 
laser-treated  films.  However,  different  regions  of  such  films  can  be  expected 
to  differ  In  microstructure  depending  on  their  rates  of  transformation  and 
therefore  on  their  thermal  history.  This  suggests  that  the  periodic 
structural  features  observed  on  laser-crystallized  films  (see  Figs.  3  and  4) 
can  be  attributed  to  oscillations  In  S  like  those  Implied  by  the  curves  for 
n  *  0.3  and  0.6  In  Fig.  2.  Furthermore,  for  sufficiently  high  values  of  n 
these  oscillations  produce  large  fluctuations  In  the  rate  of  heat  liberation 


and  therefore  In  temperature.  This  can  explain  the  periodic  fluctuations  In 
light  emission  observed  during  some  laser  crystallization  experiments  (see 
Sec.  II  C). 

B.  Comparison  with  Experiment  * 

In  order  to  carry  out  a  semi -quantitative  test  of  the  model,  we  have 
measured  the  spatial  period  of  the  structural  features  of  laser-crystallized 
Ge  films  as  a  function  of  T^. 7  Experiments  were  performed  on  amorphous  films 
0.3  urn  thick,  deposited  on  fused  silica  substrates  and  scanned  at  v  *  0.5 
cm/sec  with  a  silt  Image  of  a  cw  Nd : TAG  laser.  Initially,  a  film  at  room 
temperature  was  Irradiated  at  a  laser  power  level  just  high  enough  to  produce 
crystallization,  which  yielded  structure  In  the  transformed  film  with  a 
spatial  period  of  ~  50  gn.  In  the  following  experiments  each  film  was  heated 
to  a  successively  higher  value  of  Tt>,  the  laser  scanned  at  the  same  power 
level,  and  the  spacing  measured  after  crystallization.  This  procedure  was 
continued  until  T5  approached  the  value  resulting  In  runaway. 

The  periodic  features  obtained  by  crystallization  of  a  film  with  Tj,  of 
room  temperature  are  shown  In  Fig.  3,  an  optical  transmission  micrograph. 
These  features  are  shown  at  higher  magnification  by  the  left  side  of  Fig.  4, 
which  Is  a  brlght-fleld  micrograph  obtained  by  transmission  electron 
microscopy  (TEH)  using  125  keV  electrons.  Each  feature  consists  of  four 
different  regions:  first  a  narrow  amorphous  region,  then  a  region  containing 
a  mixture  of  amorphous  material  and  fine  grains,  next  a  broad  region  of  fine 
grains,  and  finally  another  broad  region  of  much  larger,  elongated 
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Fig.  4.  Left:  Bright-field  transmission-electron  micrograph  of  a 
laser-crystallized  Ge  film,  illustrating  four  different  microstruc¬ 
ture  regions.  Right:  Transmission-electron  diffraction  patterns 
for  fine-grained  and  large-grained  regions  (upper  and  lower,  respec¬ 
tively)  . 


crystallites  aligned  parallel  to  each  other.  The  fine-grained  region  yields 
transmission  electron  diffraction  patterns  like  the  one  shown  at  the  upper 
right  of  Fig.  4,  with  the  rings  typical  of  polycrystalline  material.  The 
large  crystallites  of  the  final  region  are  clearly  visible  as  ribbon-like 
structures  in  the  lower  left  corner  of  the  TEM  micrograph,  and  yield 
characteristic  single-crystal  transmission  electron  diffraction  patterns,  as 
Illustrated  at  the  lower  right  of  Fig.  4.  It  should  be  emphasized  that  these 
large  aligned  crystallites  are  produced  without  the  presence  of  a  relief 
structure  in  the  amorphous  silica  substrate. *2  The  nature  of  the  crystalline 
structure  within  each  period  Is  discussed  further  In  Sec.  V. 

To  use  the  model  to  calculate  the  spatial  period  In  the  film,  we  assume 
that  this  period  Is  equal  to  the  distance  AY  traversed  by  the  a-c  phase 
boundary  during  each  of  Its  successive  jumps,  from  the  point  where 
crystallization  Is  Initiated  by  the  approaching  laser  to  the  point  where  the 
boundary  comes  to  rest  ahead  of  the  laser.  This  distance  Is  just  equal  to 
(aVAx),  where  V  Is  the  normalized  laser  scanning  velocity  and  Ax  Is  the 
normalized  time  Interval  between  the  beginnings  of  two  successive  jumps,  i.e. , 
the  period  of  the  oscillations  In  S  Illustrated  by  the  curves  for  n  =  0.3  and 
0.6  In  Fig.  2. 

To  obtain  a  relationship  between  AY  and  T|> ,  we  first  used  Eq.  (3)  to 
calculate  VAx  as  a  function  of  n  for  representative  values  of  o,  V,  and  r. 
Because  of  the  Introduction  of  the  new  damping  parameter  r,  the  results 
obtained  are  numerically  somewhat  different  from  those  reported  earlier,7*8 
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but  are  qualitatively  unchanged.  The  calculated  curve  for  o*  6  n 
(corresponding  to  fixed  AT^),  V  ■  0.3,  and  r  ■  30  Is  shown  In  Fig.  5.  With 
Increasing  n.  At  and  therefore  Vat  Increase  rapidly,  leading  to  boundary 
runaway  by  n  *  0.69.  For  a  given  value  of  n.  Vat  Is  found  to  be  quite 
Insensitive  to  either  a  or  V,  showing  that  the  boundary  jump  distance  Is 
determined  primarily  by  the  properties  of  the  film  and  by  Tb,  and  does  not 
depend  strongly  on  either  the  power  or  velocity  of  the  laser. 

In  order  to  compare  theory  with  experiment,  the  curve  of  Vat  vs  n  shown 
In  Fig.  5  was  used  to  obtain  curves  relating  the  ratios  aY/aY0  and  Tb/Tc, 
where  aY0  Is  the  spatial  period  In  the  film  (l.e.,  the  boundary  jump  distance) 
for  Tb  of  room  temperature.  To  calculate  these  curves,  Eq.  (4b)  was  rewritten 
In  the  lumped-parameter  form  rpr^/(l-Tb/Tc),  where  r\j  s  fL/i4/2CTc.  Tb  and 
Tc  are  expressed  In  °C,  room  temperaure  Is  much  less  than  Tc,  so  that  no  *  n 
at  room  temperature,  and  aY0  corresponds  to  r\).  Curves  of  aY/aY0  vs  Tb/Tc  * 
1-no/n  were  calculated  by  adopting  pairs  of  numerical  values  of  m  and  Tc, 
then  compared  with  the  experimental  points  (for  each  point,  the  adopted  value 
of  Tc  was  used  to  determine  Tb/Tc).  A  reasonable  overall  fit  has  been 
obtained,  as  shown  In  Fig. 6,  for  no  *  0.22  and  Tc  -  680°C.  The  abrupt  Increase 
In  aY/aYq  when  Tb/Tc  exceeds  about  0.6  Is  associated  with  the  approach  of  Tb 
to  the  value  Tr  above  which  laser  Irradiation  results  In  boundary  runaway. 

From  Fig.  6,  using  Tc  -  680°C,  Tr  ~  [0.7  Tc]  *  500°C.  It  Is  Interesting  to 
note  that  Bagley  and  Chen10  predict  a  transformation  of  a-Ge  to  the  liquid 
state  at  a  temperature  of  696°C,  well  below  the  crystal 11ne-to-11 quid 


transition  temperature  (937°C),  and  close  to  the  value  of  ~  680°C  obtained 
from  our  analysis. 

From  the  definition  of  t^,  fL  *  wl/2cTcr*).  Taking  C  *  0.08  cal/g°C  for 
amorphous  Ge  (Ref.  3)  and  the  values  of  Tc  and  used  for  Fig.  6,  we  obtain 
fL  »  21.4  cal/g.  In  calorimetric  measurements  made  during  the  rapid  heating 
of  amorphous  Ge  films.  Fan  and  Anderson4  observed  a  sharp  transition  at  ~ 

501°C  and  measured  L  =  39.8  cal/g.  Using  this  value  for  L  yields  the  result  f 
~  0.54.  The  value  of  f  obtained  seems  large,  and  Is  discussed  further  In 
Sec.  1 1  IB. 

The  values  of  S(t)  plotted  in  Fig.  2  for  the  regions  In  which  the  phase 
boundary  moves  away  from  the  laser  Image  are  not  quantitatively  correct,  since 
the  calculated  values  of  S  In  this  region  Increase  more  steeply  as  the  size  of 
the  Interval  5t  chosen  for  the  numerical  Integration  of  Eq.  (3)  Is  decreased. 
An  examination  of  the  Integral  equation  Indicates  that  In  the  beginning  of 
these  regions  Its  solutions  Is  singular  and  the  Initial  velocity  Is  Infinite. 
This  represents  a  limitation  of  the  model  since  we  have  Implicitly  assumed 
that  growth  of  the  crystalline  phase  can  follow  at  all  times  the  requirements 
of  the  solution  of  the  heat  flow  problem.  In  fact,  the  failure  of  this 
assumption  may  account  for  the  presence  of  amorphous  and  fine-grained  material 
in  the  Initial  portion  of  each  periodic  feature  of  the  laser-crystallized 
films.  Since  the  calculated  values  of  S  at  which  the  boundary  comes  to  rest 
approach  a  limit  as  6t  decreases,  the  computed  spatial  periods  are  reasonably 


reliable. 
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The  Gaussian  width  of  the  laser  temperature  profile  can  be  determined 
from  the  relationship  a  =*  aY/VAt.  From  the  measured  AY0  »  50  nn  at  room 
temperature  and  Vat  -  0.3  corresponding  to  >  0.22  In  Fig.  5,  a  Is  ~  170  un, 
a  reasonable  value.  The  scanning  velocity  v  corresponding  to  the  normalized 
velocity  V  In  our  calculation  can  be  obtained  by  using  the  relation 
v  =  4icV/a.  A  value  for  tc  of  ~  0.09  cm2/sec  Is  estimated  by  taking 
K  =  0.035  cal/cm  sec°C,13  p  ~  5  g/cm3,  and  C  ■  0.08  cal/g°C.  The  value  of 
normalized  laser  scan  velocity  used  In  our  calculations,  which  was  chosen 
for  convenience  In  numerical  Integration,  was  V  =  0.3.  This  corresponds  to  a 
laser  scan  velocity  v  of  ~  6  cm/sec,  much  higher  than  the  actual  value  of  0.5 
cm/sec  used  in  our  experiments.  As  indicated  above,  the  results  of  our 
analysis  of  periodic  a-c  boundary  motion  are  Insensitive  to  the  value  of  v 
assumed  for  the  calculation  as  long  as  v  Is  much  less  than  the  boundary 
velocity  vac.  For  both  runaway6  and  periodic  motion9  of  the  boundary, 
experimental  values  of  vac  lie  In  the  range  of  100  -300  cm/sec. 

C.  Temperature  Distribution  During  Periodic  Boundary  Motion 

Once  the  a-c  boundary  position  as  a  function  of  time,  S(t),  Is  known, 
Eq.  (1)  can  be  used  to  calculate  the  temperature  at  an  arbitrary  point  S" 

In  the  film  at  time  t.  The  resulting  expression  for  the  reduced  temperature 
In  terms  of  normalized  quantities  In  the  frame  of  reference  moving  with  the 
laser  Is: 
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T'  (S',x)  *  a  exp  {-[S']2} 


+  n  /T  C$  (t')  +  V]/(t-t')1/2>  x 
o 

exp  {-r  (T-x')}j|exp  {-[S'  -  S(t')  +  V 
where  the  reduced  temperature  Is 
T'  (S'.t)  =  T(S',t)-  Tb 

fc  -Tb" 


dx' 


(5) 


It  can  be  seen  that  when  S'  Is  the  position  of  the  a-c  boundary,  S(x), 

T(S',x)  =  Tc,  T'(S',x)  =  1,  and  Eq.  (5)  becomes  the  Integral  equation  (3). 

The  numerical  evaluation  of  Eq.  (5)  Is  discussed  In  Appendix  A. 

To  illustrate  the  temperature  fluctuations  occurring  during  periodic 
boundary  motion,  we  have  used  Eq.  (5)  to  calculate  T'(S', x)  as  a  function  of 
S'  for  several  fixed  values  of  x,  for  the  case  of  3.6,  n*  0.6,  V  =  0.3,  and  r 
*  30.  The  results  are  shown  In  Fig.  7.  It  can  be  seen  that  during  the  forward 
motion  of  the  a-c  phase  boundary,  a  temperaure  pulse  develops  and  propagates 
away  from  the  center  of  the  laser  temperature  profile.  Before  the  pulse  can 
escape,  however.  Its  motion  stalls  and  It  decays. 

The  case  of  the  periodic  fluctuations  In  light  emission  observed  In 
some  experiments  Is  Illustrated  In  Fig.  8  by  the  case  crl.2,  r\  *  0.6,  V  *  0.3, 
and  r  *  30.  In  this  case  the  laser  temperature  contribution  Is  just  barely 
large  enough  to  Initiate  a-c  boundary  motion,  and  the  latent  heat  released  Is 
enough  to  cause  a  sharp  rise  In  temperature  accompanied  by  the  emission  of 
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black  body  visible  radiation.  The  peak  temperature  Tp  reached  can  be 
calculated  for  this  example  by  noting  that 

Tp  -  Tb  +  TP'  (Tc  -  Tb),  (6) 

where  Tp'  Is  the  peak  value  reached  In  Fig.  8.  Equation  (6)  can  be  rewritten 

Tp  -  Tc  [(1-  ho)  +  Tp'(J)o)]  . 
n  n 

Using  Tc  »  680°C,  *  0.22,  n  *  0.6,  and  Tp'  -  2.2,  we  have 

Tp  -  986°C  . 

Although  the  numerical  result  Is  only  approximate  (since  the  computed 
velocity  of  the  a-c  boundary  depends  on  the  time  Interval  Sr  chosen  for 
numerical  Integration),  the  calculation  Indicates  that  a  large  rise  In 
temperature,  and  therefore  light  emission,  can  occur  during  a-c  boundary 
motion. 

III.  Steady  State  Runaway  Solutions 
A.  Integral  Equation  Model 

As  we  have  seen  In  Sec.  IIA,  for  large  enough  values  of  n  and 
relatively  small  values  of  normalized  laser  scan  velocity  V,  the  solutions  of 
Eq.  (3)  are  of  the  runaway  type  with  the  a-c  boundary  escaping  from  the 
region  of  the  laser  Image  and  moving  far  ahead  of  It.  In  this  case,  Eq.  (3) 
can  be  solved  exactly  without  the  need  for  numerical  Integration,  by  making 
the  assumption  that  a  steady  state,  constant  velocity  solution  exists  after 
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the  passage  of  a  long  enough  time  tfc.  For  the  boundary  far  ahead  of  the  laser 
Image  the  laser  contribution  to  Eq.  (3)  Is  negligible.  We  rewrite  Eq.  (3) 
at  time  t0  In  the  frame  of  reference  moving  with  the  boundary  at  the 
normalized  velocity  Vac  »  V,  with  Vac  to  be  determined.  In  this  frame  of 
reference,  V  Is  replaced  by  Vac,  S(t)  -  S(t'),  and  §(T')  -  0  (since  the 
solution  Is  stationary).  Equation  (3)  then  becomes 


1  =  n  J  [Vac/(T-  t')1/2]  x 


')j^xp  {-[Vac  (t-t')]2/(t-t')}  dT'.  (7) 


jxp  {-r  (t-t 


Introducing  the  new  variable  z2  =  (r  +  Vac2)  (t-t'),  Eq.  (7)  can  be  rewritten 


1  *  n 


[r+  Vac2] 


,  [r+  Vac2]l/2Tol/2  _z2 

x  2  /  e  dz 

o 


n  «l/2  vac  erf  {[r  +  Vac2]l/2^1/2  j 
[r  ♦  vac2]  72 


where  erf  {x}  »  2  yx  e-z  dZ>  error  function. 

/i  o 
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The  error  function  approaches  1  for  x  >  2.  Therefore,  for 
[r  +  V^2]1/2  '*b  >  2»  Ed*  (8)  reduces  to 

1  =  n w1/2  vac  .  (9) 

- TT75" 

[r+  vac2] 

Equation  (9)  can  be  used  to  determine  the  runaway  a-c  boundary  velocity  Vac 
from  the  values  of  n  and  r: 

Var  *  r1/2  .  (10) 

T72 

[  n?  n-1] 

It  can  be  seen  that  Eq.  (10)  has  no  solution  for  any  value  of  Vac,  no  matter 
how  large,  unless  n  w  1/2  >  1,  or  n  >  0.56, 

In  the  runaway  case,  Eq.  (5)  for  the  reduced  temperature  distribution 
can  also  be  evaluated  exactly  for  large  iq.  In  the  same  frame  of  reference 
moving  with  normalized  velocity  Vac,  Eq.  (5)  becomes 

T'(S')  -  n  JT°  [Vac/( t-t')1/2] 
o 

xp  {-[S'+  Vac  (t-t')]2  /  (T-T')}|dT'  ,  (11) 
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xp  {-r(T-T')> 


where  S'  is  measured  from  the  position  of  the  a-c  phase  boundary.  Eg.  (11) 
can  be  written. 


T'(S')  = 


[r  ♦  vac2] 


[r  +  Vgc2]1/2  t01  n 


2  e-2S 


Vac/° 


2  2,2 

e-z  e-  Q  H  d2, 


where  q2  =  S'2  [r  +  Vac2]  . 

Using  the  known  definite  integral, 


.  _p2z2  _q2/z2  -2pq 

/  e  e  dz  =  /tr  e  ,  (p,q  >  0) 

o  Tq 

we  obtain,  for  large  x0. 


T'(S')»  nitl/2  i/g  Vac 
[r+  vac2] 


e  -2S'Vac  e-2  Cr  +  Vac2]l/2  IS'1  .  (13) 


For  S'  s  0,  corresponding  to  the  position  of  the  a-c  boundary,  T'(S')  =  1,  and 
Eq.  (13)  reduces  to  Eq.  (9),  as  It  should. 

To  explore  the  approach  to  the  steady-state  runaway  temperature 
distribution,  we  have  used  the  values  of  S(t)  obtained  from  numerical 
Integration  of  Eq.  (3)  for  a  ■  2.0,  n  *  0.8,  V  ■  0.3,  and  r  •  30,  to  calculate 


27 


T'(S',t)  from  Eq.  (5)  for  several  values  of  t.  The  results  are  shown  In  Fig. 
9.  In  contrast  to  the  periodic  solution  case  shown  In  Fig.  7,  the  a-c 
boundary  in  the  runaway  case  develops  a  well-defined  peak-structure  In 
T'(S',t),  which  propagates  away  from  the  laser  Image  as  r  Increases.  The 
velocity  of  propagation  of  the  a-c  boundary  and  the  amplitude  and  shape  of  the 
propagating  temperature  pulse  shown  In  Fig.  9  are  not  In  agreement  with  the 
exact  limiting  results  given  by  Eqs.  (10)  and  (13),  respectively.  Thus,  from 
Eq.  (10)  with  r  =>  30  and  n  3  0.8  we  have  the  exact  normalized  velocity  Vac  * 
5.45,  much  less  than  the  normalized  velocity  obtained  from  numerical 
Integration  and  Indicated  by  Fig.  9.  These  discrepancies  may  In  part  be  due 
to  computational  problems,  but  may  also  be  due  to  a  slow  approach  to  the 
steady  state.  (It  should  be  noted  that  at  the  maximum  S(t)  shown,  the  a-c 
boundary  has  moved  only  a  few  times  the  Gaussian  width,  ~  170  tin,  from  the 
laser  beam.) 

B.  Two-Dimensional  Solution 

In  Appendix  B  we  obtain  expressions  for  the  temperature  distribution  In 
the  steady  state  a-c  boundary  runaway  condition  for  a  semiconductor  film  as  a 
function  of  normalized  distance  S'  from  the  boundary  by  the  solution  of  a 
two-dimensional  bounary  value  problem.  In  this  treatment  the  flow  of  heat 
Into  the  substrate  Is  handled  exactly,  and  It  Is  not  necessary  to  Introduce 
ad-hoc  parameters  equivalent  to  f  and  y  of  the  one-dimensional  Integral 
equation  model.  The  solution  of  the  problem  Is  discussed  In  Appendix  B. 


Equations  (B31a)  and  (B31b)  give  the  temperature  distribution  Ti(S')  In 
the  semiconductor  film  of  thickness  bi  supported  on  a  substrate  of  thickness 
b2  as  a  function  of  steady  state  velocity  V  of  the  a-c  boundary.  If  the  a-c 
boundary  Is  at  a  temperature  Tc,  then  we  can  rewrite  Eqs.  (B31a)  and  (B31b)  In 
the  form 


...  r  -2VacS'  -2[rn+  +  Vac2]l/2  |S'| 

T'(S')  *  n"Vac  wl/2  l  e  e _  _  f  (Vac)  ,(S'  <  0) 


n+ 


[  rn+  ♦  Vac2] 


T7T 


(14a) 


„  -2VacS '  -2[rn_  +  Vac2]l/2  |S'| 

T "(S ")  *  n'Vjc  nl/2  i  e  e _  _  f  (Vac),  (S'  >  0)  (14b) 

n-  „  1/2  n- 

[rn-  ♦  Vac2] 

where 


T'(S')  *  1^1  (S')  -  Tb  , 
"  Tc~  Tb  ’ 


L _  .  (14c) 

T/2 

*  C!(Tc  -  Tb) 


f  (V) 

n± 


1-6/2  (cot  Qn±  .  csczQn±l 

(l  ±  V(M-l)  ) 

1  Qn*  ' 

t  /(rn±+vz)i. 

B  *  K2b2/*lbl  and  M  »  q/t?  ,  and  S'  Is  the  normalized  distance  from 
the  a-c  boundary,  Vac  Is  the  normalized  a-c  boundary  velocity.  For  the 
normalization  unit  of  length  d  we  use  b2.  Subscripts  1  and  2  refer  to  laser 
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crystallized  film  and  substrate,  respectively.  The  rn+  and  rn_are  roots  of  a 
transcendental  equation,  (826b),  and  the  Qn±  are  given  by  the  expression 
(B26a)  evaluated  for  rn±. 

One  case  for  which  the  solution  of  Eqs.  (14a)  and  (14b)  can  be  readily 
evaluated  1$  a  thin  amorphous  film  on  a  thick  substrate  with  Identical  thermal 
properties  (M  =  1,  $  *  fy/bi  »  1).  For  this  case, 

f  =  f  =  2/0;  r  ■  r  =  (2n  +  1)2,2/16. 
n+  n-  n+  n" 


The  general  condition  determining  the  steady  state  a-c  boundary 
velocity  Vac  Is  that  at  S'  =  0,  T'(0)  =  1  or. 


1  *  n'Vac,l/2  £ 
n+ 


n  +  (Vac) 

"  1/2 

[r  +  Vac2] 


n'Vac'ir1/2  l 
n- 


n-  (Vac) 

,  1/2 

[r  +vac2] 


(15) 


It  Is  Interesting  to  compare  the  results  of  Eqs.  (14)  and  (15)  with  the 
one-dimensional  Integral  equation  results  of  Eqs.  (13)  and  (9),  respectively. 
Noting  that  n  ■  fn'  ,  It  can  be  seen  that  EqSo  (13)  and  (9)  represent 
approximations  In  which  single  average  terms  have  replaced  sums  of  terms  of 
the  same  form  In  the  two  dimensional  result. 
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We  define  the  functions  In  the  summation  In  Eq.  (14a)  and  (14b)  as 


F^S'.V)  =  l 
n  ± 


-2VS'  -2[r  +  v2]l/2  |S'| 

e  e  n± 

[r  +  V2] 
n* 


f  (V). 

n± 


(16) 


F±(S',V)  Is  proportional  to  the  steady-state  temperature  distribution  due  to  a 
slit  Image  source  of  heat  scanned  with  normalized  velocity  V  along  the  surface 
of  a  semiconductor  film.  When  evaluated  for  V  =  Vac,  It  yields  the  shape  of 
the  distribution  for  a-c  boundary  runaway.  We  have  evaluated  Eq.  (16)  as  a 
function  of  S'  for  a  number  of  values  of  V  for  the  case  of  a  thin  film  of  Ge  on 
a  fused  silica  substrate,  since  the  experimental  results  reported  in  Sec.  II  B 
were  obtained  for  this  case.  The  values  of  Kj  =  0.035  cal /cm  sec°C,  icj  *  0.09 
cm2/sec  and  K2  =  0.0025  cal /cm  sec°C,  *2  “  0.005  cm2/sec  were  those  for  Ge 
(1),  and  fused  silica  (2).  The  thicknesses  used  in  the  calculation, 
appropriate  to  our  experiments7  were  bj  =0.3  vm  for  Ge  and  b2  =  1mm  for  fused 
silica.  With  these  values,  6  2  240,  and  M  2  18. 

Figure  10  shows  F(S',V)  as  a  function  of  S'  for  several  values  of  V. 

For  V  J'O.l,  the  leading  edge  of  the  temperature  distribution  begins  to 
sharpen,  and  the  peak  value  begins  to  drop.  This  Is  an  Indication  that  V  has 
reached  a  value  at  which  diffusion  of  heat  can  no  longer  restore  the  static 
temperature  distribution  as  the  heat  source  moves.  It  should  also  be  noted 
that  the  fact  the  two  portions  of  the  curves  of  F*  (S',V)  meet  at  S'  =  0  is  a 
significant  check  on  the  correctness  of  the  numerical  calculation,  since  the 


Fig.  10.  F($',  v)  as  a  function  of  S'  for  several  values  of  V. 
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formal  expressions  for  [F+(S',V)]s-  >  o  and  [F*(S',V)]$'  .  o  are  not 
Identical . 

To  obtain  quantitative  results  for  Vac,  we  make  use  of  Eq.  (15)  and 
assume  that  Tc  Is  a  fixed  transition  temperature  (modification  of  this 
assumption  Is  discussed  In  Sec.  V).  The  relation  to  be  solved  for  V  can  be 
written, 

x1/2n'Vac  -  [F(0,  Vac)]-1  .  (17) 

We  plot  log  [F(0,V)]-1  vs  log  V  and  log  [w1/2n'V]  vs  log  V  on  the  same  graph. 

The  solution  for  Vac  will  be  given  by  the  Intersection  of  the  two  curves. 

Figure  11  shows  [F(O.V)]-1  for  H  =  18,  and  for  several  values  of  8.  We 
consider  first  the  curve  for  6  *  240,  corresponding  to  our  experiment.  One  of 
the  straight  lines  shown  gives  log  [i4/2n'V]  for  n'  ■  2.2,  and  we  see  that 
there  Is  no  solution  for  this  value  of  n".  The  minimum  value  of  n'  for  a 
steady  state  solution,  n'mlnt  corresponds  to  the  lowest  temperature  (Tr)  at 
which  runaway  can  occur.  This  solution  Is  obtained  when  the  two  curves  first 
cross,  and  Is  found  to  occur  when  n'  a  6.1.  For  values  of  r\'  greater  than 

•% 

6.1,  there  are  steady  state  or  runaway  solutions  for  Vac  which  decrease  with  $ 

* 

Increasing  n".  It  should  be  noted  that  the  two  curves  first  make  contact  I 

tangentially  as  n'  Increases  to  n'min«  This  Is  similar  to  the  Integral  ‘  ? 

equation  result,  Eq.  (10).  j 

It  Is  difficult  to  Identify  the  first  true  "crossing"  point  for  two  j 

I 

curves  that  approach  each  other  tangentially.  Roughly,  for  n'  a  6.1  the  two  i 
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Fig.  11.  Plots  of  log  [F(0,  V)]  vs  log  V  for  8  -  240,  120,  and  60 
plots  of  log  [ir1/Z  n'V]  vs  log  V  for  n'  =  9.0  and  tT  =  2.2. 
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curves  cease  to  run  parallel  for  log  V  ~  1.6.  Using  the  relation  v  ~  4ieiV/b2, 
with  k\  *  0.09  cm2/sec,  and  b2  =  10"1  cm,  we  find,  vac  s  140  cm/sec.  Larger 
values  of  n'  should  lead  to  smaller  values  of  vac,  consistent  with  the 
integral  equation  resUlt,  Eq.  (10). 

From  the  relation  (14c),  we  can  calculate  the  value  of  Tc  If  we  equate 
t)'  with  o'mln  =  6*1  and  Tb  with  the  experimental  value  of  Tr,  ~  500°C.  This 
leads  to  Tc  ~  550°C,  much  lower  than  the  value  of  ~  680°C  obtained  from  the  fit 
of  the  periodicity  data  to  the  theoretical  results  obtained  from  the  Integral 
equation  model.  As  we  shall  see  shortly,  the  discrepancy  can  be  explained 
qualitatively  by  the  presence  of  a  thermal  barrier  between  film  and  substrate. 

We  can  attempt  to  approximate  the  calculated  F(0,V)  In  the  form 

F (0,  V)  a  f/[r  +  v2]l/2  ,  (17a) 

where  f  and  r  are  to  be  determined  for  the  best  fit.  Equation  (17)  would  then 
be  precisely  of  the  form  Eq.  (9),  obtained  from  the  Integral  equation,  where 
n=fn'.  The  open  circles  in  Fig.  11,  representing  Eq.  (17a)  with  r  *  7.95  and 
f  =  0.093,  fit  the  calculated  curve  of  log  [F(0,V)]-1  with  0  =  240  and  M  »  18 
for  small  V  and  large  V,  but  fit  poorly  In  between.  This  value  of  f  Is  much 
lower  than  that  obtained  from  the  Integral  equation  fit  to  the  periodicity 
data,  f  ~  0.58.  As  we  shall  see  shortly,  this  descrepancy  can  also  be 
qualitatively  explained  by  the  presence  of  a  thermal  barrier  between  film  and 


substrate 


We  have  examined  the  effect  of  changing  bj  and  on  the  values  of  a 
number  of  properties  of  a-c  boundary  runaway.  Figure  11  shows  plots  of  log 
CF(0,V}]_1  vs  log  V  for  M  *  18  and  8  ■  120  and  60,  In  addition  to  the  curve 
for  8  =  240,  which  corresponds  to  our  experiment.  We  can  Interpret  the 
reduction  In  0  as  due  to  either  a  decrease  In  b2  or  an  Increase  In  bj.  The 
results  obtained  for  Vac  and  vac  for  n'  ■  9.0  (corresponding  to  fixed  Tj>,  and 
well  above  n'mln  for  all  3  values  of  0)  are  listed  In  Table  I.  It  can  be  seen 
that  for  bi  held  fixed  at  0.3  un  ,  vac  decreases  only  slightly  as  b2 
decreases.  On  the  other  hand,  for  b2  held  fixed  at  103  |/n,  vac  decreases 
almost  Inversely  with  Increasing  bj.  However,  In  most  cases  It  might  be 
difficult  to  see  this  striking  drop  In  vac,  since  at  temperatures  above  Tr, 
spontaneous  fluctuations  can  cause  the  a>c  transformation  (see  Sec.  V.) 

Also  shown  In  Table  I  are  values  of  n'mln.  with  corresponding  values 
for  Vac  and  vac,  as  well  as  approximate  values  of  f  and  r.  It  can  be  seen  that 
f  Increases  nearly  Inversely  as  0  decreases,  while  r  Increases  slowly.  It 
should  also  be  noted  that  n'mln  decreases  as  bi  Increases.  This  means  that 
the  lowest  value  of  Tb  for  runaway,  Tr,  decreases  as  film  thickness  Increases. 
As  the  film  thickness  bj  Increases,  the  velocity  vac  for  n'mln  Increases,  but 
remains  In  the  range  of  100-300  cm/sec. 

The  small  values  obtained  for  Tc  and  f  by  applying  the  present  theory 
directly  to  experiment  can  be  understood  If  we  postulate  the  presence  of  a 
thermal  barrier  between  the  fused  silica  substrate  and  the  amorphous  Ge  film 
In  our  laser  crystallization  experiments.7  The  presence  of  a  thermal  barrier 
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TABLE  I 


Results  of  Two  Dimensional  Solution  for  a-c  Boundary  Notion 

Values  of  a  number  of  quantities  calculated  from  exact  solutions  of  a  model  of  laser 
crystallization,  as  obtained  from  Fig.  11.  For  three  values  of  the  parameter  0, 
effective  parameters  f  and  r  and  the  minimum  value  of  n"  for  runaway,  n'  *in,  are  given. 
With  n'  *  n"  min  and  n'  ■  9.0,  normalized  velocities  Vac  for  a-c  boundary  runaway  are 
obtained.  For  combinations  of  substrate  (62)  and  film  (bi)  thicknesses  given,  these 
values  of  Vac  lead  to  the  a-c  boundary  runaway  velocities  vac  that  are  shown. 


6 

f 

r 

n'mln 

Vac  Var 

(n'-n'min)  (n'-9.0) 

t>2 

(U») 

bl 

(um) 

vac(c*/*«c) 

(i'“»>'m1n) 

Vac(c*/*»c) 

(n'-9.0) 

240 

0.093 

7.95 

6.1 

39  10 

103 

0.3 

140 

36 

120 

0.178 

8.72 

3.2 

47  4.07 

103 

0.6 

170 

14.7 

5x10* 

0.3 

29.3 

60 

0.331 

10.15 

1.7 

71  1.9 

103 

1.2 

255 

6.8 

2.5x102 

0.3 

27.4 

Is  not  at  all  unreasonable,  and  could  be  readily  Included  In  the  derivation  of 
modified  equations  (14a)  and  (14b).  The  effect  of  a  thermal  barrier  Is 
approximately  equivalent  to  decreasing  the  thermal  conductivity  of  the 
substrate,  and  Is  therefore  qualitatively  described  by  decreasing  0  (although 
the  thermal  conductivity  ratio  appears  In  M,  this  describes  the  way  heat 
diffuses  once  It  Is  transferred  to  the  substrate  and  should  therefore  be  kept 
fixed  for  this  crude  argument).  To  obtain  a  rough  Idea  of  the  effect  of  such 
a  barrier,  we  examine  the  effect  of  decreasing  0  by  comparing  the  results  for 
0  =  240  with  those  obtained  for  0  =  60.  Referring  to  table  I,  we  see  that  for 
0  =  60,  f  has  Increased  to  0.331,  much  closer  to  the  empirical  value  0.58 
determined  from  the  experimental  data  by  comparison  with  the  one-dimensional 
Integral  solution.  Furthermore,  n'  min  has  dropped  from  6.1  to  1.7,  yielding 
a  value  of  Tc  s  665°  determined  from  equation  (14c).  This  Is  close  to  the 
value  of  680°  C  obtained  from  the  fit  of  experimental  data  to  the  Integral 
equation  model.  While  this  argument  Is  very  crude.  It  does  suggest  that  the 
presence  of  a  thermal  barrier  can  eliminate  the  discrepancy  between 
experimentally  determined  parameters  and  the  results  of  a  first-principles 
calculation. 

IV.  Steady  State  Laser-Guided  Solutions 

An  Intriguing  class  of  steady  state  solutions  exists  In  which  the  a-c 
boundary  proceeds  at  a  constant  velocity  on  the  leading  edge  of  the  scanning 
laser  Image.  Solutions  of  this  type  will  be  obtained  for  both  the  Integral 
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equation  model  and  the  two-dimensional  model  of  the  a-c  boundary  motion.  The 
laser-guided  mode  of  a-c  transformation  could  have  Important  Implications  for 
the  controlled  growth  of  semiconductor  films.  The  question  of  the  Initial 
conditions  which  permit  the  system  to  settle  Into  the  steady  state 
laser-guided  mode  of  motion  will  be  discussed  In  Sec.  V. 

A.  Integral  Equation  Model 

To  find  laser-guided  solutions,  we  rewrite  Eq.  (3)  In 
the  frame  of  reference,  moving  at  normalized  velocity  V,  In  which  the  laser 
and  a-c  boundary  temperature  distributions  are  stationary.  The  result  Is 


1  *  T^(0)  +  n  /  [V/fT-r')1/2]  [exp  {-r 

x[exp  {-[V  (T-T')]2/(T-T'))]dT'  , 


(18) 


where  TA'(0)  Is  the  reduced  temperature  due  to  the  scanning  laser  Image  at 
the  poslton  of  the  a-c  boundary,  T4'(0)  *  T*(0)/(Tc-Tb). 

As  In  Sec.  Ill  A,  the  Integral  In  Eq.  (18)  can  be  evaluated  exactly  for 
large  enough  -iq,  and  we  find 

1  *  Tt'(0)  +  hVw 1/2  /  [r  +  v2]l/2  .  (19) 

According  to  Eq.  (19),  for  values  of  n  and  V  such  that  r^irl/2/[r  +  v2]l/2  <  i, 
there  Is  a  value  of  T ^(O)  which  will  satisfy  the  condition  for  a 
laser-guided  solution.  For  given  values  of  V  (determined  by  laser  scan 
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velocity)  and  n  and  a  given  reduced  temperature  profile  due  to  the  laser,  the 
(a-c)  boundary  will  ride  at  a  distance  In  front  of  the  laser  Image  such  that 
Eq.  (19)  Is  satisfied.  It  should  be  noted  that  solutions  of  Eq.  (19)  exist 
for  values  of  V  arbitrarily  small.  The  question  of  coupling  Into  such 
solutions  In  real  situations  Is  discussed  In  Sec.  V. 

B.  Two-Dimensional  Solution 

It  Is  a  straightforward  matter  to  show  that  the  two-dimensional 
solution  of  the  laser  crystallization  problem  also  leads  to  laser-guided 
stationary  solutions.  As  shown  In  Appendix  B,  the  solution  of  the  steady 
state  temperature  distribution  for  the  boundary  value  problem  In  the  presence 
of  several  heat  sources  Is  the  superposition  of  the  solutions  In  the  presence 
of  the  Individual  heat  sources.  In  the  presence  of  both  a  laser  source  of 
heat  and  the  heat  source  due  to  the  moving  a-c  boundary,  the  solution  for 
steady  state  runaway,  Eq.  (15),  becomes 


1 


T*'(0)  +  n^V*1/2  I 
n+ 


fn+  (V) 

—  rw 

[rn+  +  V2] 


-  Tt'(0)  +  n-Vxl/2  l 
n- 


nd^L 

rn-+  V2] 


(20) 


where  TA'(0)  Is  the  reduced  temperature  due  to  the  scanning  laser  Image  at  the 
position  of  the  a-c  boundary.  The  close  correspondence  between  the 
two-dimensional  solution  Eq.  (20)  and  the  Integral  equation  solution  Eq.  (19) 
should  be  noted. 
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Uslng  the  definition  of  F*  (S',V),  Eq.  (16),  we  obtain  for  the 
laser-guided  problem  an  expression  analogous  to  Eq.  (17),  which  was  solved 
to  determine  Vac  for  runaway.  The  result  Is 


wl/2n-V  -  [F(0,V)]-1  .  (21) 

i-Mo) 

We  may  think  of  n7 [l-T a"(0) ]  as  an  effective  value,  n'eff*  Referring  to  Sec. 
Ill  B  and  Fig.  11,  It  can  then  be  seen  that  by  decreasing  l-T^'tO),  the 
normalized  laser  scan  velocity  V  for  laser-guided  steady  state  solutions  can 
be  reduced  greatly.  For  a  given  reduced  temperature  distribution  due  to  the 
laser  scanning  at  velocity  V,  the  a-c  boundary  rides  at  a  point  In  front  of 
the  laser  Image  satisfying  Eq.  (21). 

V.  Discussion 

One  of  the  Important  assumptions  made  In  calculating  the  motion  of  the 
a-c  boundary  for  both  the  one-dimensional  Integral  equation  and 
two-dimensional  solution  was  the  existence  of  a  sharply  defined  transition 
temperature  Tc.  In  fact,  the  velocity  of  growth  of  crystalline  material 
directly  from  the  amorphous  state,  §(T),  Is  a  function  of  temperature. 
Increasing  rapidly  as  the  temperature  Is  raised.  This  could  be  Interpreted  as 
determining  the  temperature  TC(S)  at  which  the  transition  can  take  place  at  a 
value  of  $  demanded  by  the  heat  flow  calculation.  From  this  simplistic  point 
of  view,  over  a  wide  range  of  values  of  §,  Tc  varies  relatively  little  and  can 
be  taken  as  constant.  As  already  pointed  out  In  Sec.  II  B,  the  model  solution 
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of  the  Integral  equation  at  the  beginning  of  each  period  of  boundary  motion 
requires  Infinite  velocity,  and  the  presence  of  amorphous  and  fine-grained 
material  In  the  initial  portion  of  each  period  may  be  related  to  the  limits  on 
the  rate  at  which  crystalline  growth  can  occur. 

We  could  attempt  to  Improve  the  approximate  description  of  the  boundary 
motion  using  the  Integral  Eq.  (3)  by  replacing  n  by  the  quantity 

neff(S)  =  f  L(S)/[TC(S)  -Tb]  .  (22) 

where  TC(S)  is  the  critical  temperature  appropriate  to  the  value  of  S,  and 
L(§)  is  an  effective  latent  heat  which  Is  that  portion  emitted  when  partial 
conversion  to  crystalline  material  takes  place  at  boundary  velocity  $.  We  now 
Interpret  §  as  the  velocity  of  a  "heat  flow  boundary",  which  Is  identical  to 
the  a-c  boundary  when  that  boundary  Is  well  defined.  We  can  see  that  at  large 
§,  ieff(S)  Is  reduced  because  L(§)  Is  reduced  and  TC(S)  Is  Increased.  At 
lower  values  of  §,  L(§)  ♦  L,  and  TC(S)  ~  Tc.  The  Introduction  of  an  neff(S) 
of  this  form  could  eliminate  the  unphysical  singularity  In  the  Initial  motion 
present  In  Eq.  (3),  while  still  yielding  a  very  high  Initial  value  of  §, 
consistent  with  the  Initial  amorphous  and  fine-grained  polycrystalline 
material  In  each  period.  The  Introduction  of  an  n'eff(vac)  E<1*  U7)  would 
also  modify  the  exact  solutions  obtained  for  Vac  from  the  boundary  value 
problem. 

As  mentioned  In  Sec.  II  B,  the  value  of  Tc  obtained  by  fitting  the 
results  of  the  Integral  equation  solution  to  our  experimental  data  was 
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Tc  -  680°C,  very  close  to  the  value  of  696°C  predicted  by  Bagley  and  Chen10 
for  the  transition  temperature  from  the  amorphous  to  liquid  state  of  Ge.  It 
has  been  suggested14*15  that  "explosive"  crystallization  (or  runaway)  may 
correspond  to  the  occurrence  of  this  transition,  and  some  experimental 
evidence  supporting  this  point  of  view  has  been  obtained.15  If,  in  fact,  Tc 
does  correspond  to  the  amorphous-to-1 iquld  transition,  the  assumption  of  a 
fixed  Tc  over  a  range  of  S  Is  reasonable.  Once  the  transition  from  the 
amorphous  to  the  liquid  state  had  occurred,  the  unstable  liquid  would  rapidly 
transform  to  the  crystalline  state.10  For  small  temperature  differences 
across  a  thin  liquid  layer,  the  treatment  we  have  presented  would  then  require 
no  major  modification. 

Throughout  this  paper,  we  have  treated  the  amorphous  and  crystalline 
regions  of  the  semiconductor  film  undergoing  laser  crystallization  as  Istropic 
and  homogeneous  regions,  with  the  transformation  from  the  amorphous  to  the 
crystalline  state  governed  by  macroscopic  heat  flow  equations.  An  Important 
ingredient  absent  from  this  description,  and  essential  for  even  a  qualitative 
understanding  of  crystallized  film  morphology.  Is  the  role  of  nucleatlon 
centers17  and  "nucleatlon  events".  Because  the  exact  nature  of  the  a-c 
transformation  In  laser  crystal  1 Izatlon  Is  not  yet  certain,  we  will  not  vj 

attempt  a  detailed  description  of  nucleatlon.  We  regard  a  nucleatlon  event  as  \ 

a  spontaneous,  localized  fluctuation  which  initiates  a  transformation  from  the  _  j 

amorphous  to  the  crystalline  state  and  produces  a  small  Increase  In  local  j 

temperature  above  the  background  determined  by  macroscopic  heat  flow.  A  * 

i 

! 
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nucleatlon  event  has  a  probability  of  occurrence  which  Increases  with  the 
density  of  nucleatlon  centers  (which  depends  on  film  perfection)  and  with 
Increasing  temperature.  For  T^.  just  below  Tr,  the  lowest  temperature  for 
runaway,  nucleatlon  events  can  produce  transformation  from  the  amorphous  to 
crystalline  state  in  small  regions  of  finite  range  about  nucleation  centers.18 
For  Tb  ~  Tr  and  in  the  absence  of  an  external  disturbance  such  as  a  laser 
beam,  multiple  nucleatlon  events  can  occur,  overlapping  and  spontaneously 
transforming  the  entire  film  to  the  crystalline  state  in  a  short  time  If  the 
probability  of  such  events  Is  high  enough.  From  the  discussion  of  Sec.  Ill  B, 
the  lowest  value  of  Tb  for  runaway,  Tr,  is  higher  for  thinner  films, 
decreasing  as  the  films  become  thicker.  Since  the  probability  of  nucleatlon 
events  Increases  with  Increasing  temperature,  spontaneous  transformation  to 
the  crystalline  state  should  occur  most  readily  at  Tb  ~  Tr  for  thin  amorphous 
films.  This  Is  consistent  with  our  observation  of  laser-induced  runaway  for 
films  of  thickness  greater  than  1  yn,  while  only  spontaneous  transformation  is 
observed  for  films  0.3  un  thick. 

The  occurrence  of  nucleatlon  events  appears  to  play  a  major  role  in  the 
morphology  of  laser-crystallized  films.  In  the  films  that  we  have  examined  in 
detail,  the  elongated  crystallites  within  each  periodic  feature  form  a  roughly 
chevron-like  pattern,  with  the  two  halves  of  the  pattern  symmetrical  about  an 
axis  that  Is  parallel  to  the  laser  scan  direction  and  located  near  the  center 
of  the  laser  Image.  The  elongated  crystallites  on  each  side  of  the  pattern 
have  their  long  axes  aligned  along  [100]  directions.  A  blown-up  picture  of  a 
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crystallized  film  near  the  center  axis  of  the  chevron  pattern  Is  shown  In  Fig. 
12.  The  ''origin"  of  the  chevron  pattern  within  each  periodic  feature  is 
a  small  region  located  somewhere  near  the  center  axis.  The  exact  lateral 
position  of  the  origin  varies  from  one  period  to  the  next  over  a  distance  of 
perhaps  100  un  (~  10%  of  the  slit  Image  length).  This  suggests  that  the 
origin  represents  the  position  of  the  first  nucleation  event  within  a 
periodic  feature  where  the  heat  flow  boundary  is  moving  slowly  enough  for 
large  crystallites  to  form.  The  first  nucleation  event  would  be  likely  to 
occur  near  the  center  axis,  since  the  laser  slit  Image  is  an  ellipse  of  high 
aspect  ratio  and  the  temperature  should  therefore  first  reach  Tc  near  this 
axis.  Once  nucleation  and  growth  of  large  crystallites  began,  heat  would  flow 
forward,  but  also  laterally,  raising  the  temperature  and  Inciting  the 
nucleation  of  new  growth  centers  laterally.  The  direction  of  growth  from 
these  centers  would  tend  to  be  channeled  by  such  factors  as  Interaction  with 
neighboring  crystallites  already  formed, 19  strong  anisotropy  In  the 
directional  dependence  of  growth  rate,  and  perhaps  stresses  In  the  film.  The 
result  of  all  of  these  influences  would  combine  in  a  statistical  way,  but  it 
can  be  argued  that  the  outcome  would  be  largely  determined  by  the  tendency  of 
growth  to  occur  along  the  direction  of  the  temperature  gradient  at  the  a-c 
boundary,  combined  with  the  high  rate  of  growth  along  a  [100]  direction 
compared  to  other  crystallographic  directions.  Motion  of  the  a-c  boundary 
would  occur  laterally  as  well  as  along  the  laser  scan  direction,  but  a 
periodic  morphological  pattern  would  still  result.  This  description  Is 
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LASER  SCAN  DIRECTION  - 


- ►  100  pm  «< - 


Fig.  12.  Blown  up  view  of  periodic  features  in  a  Ge  film  crystallized 
by  the  slit  image  of  a  laser.  Picture  taken  near  the  center  axis  of 
the  slit  image  shows  the  origins  of  chevron  pattern. 


consistent  with  the  observed  chevron  pattern,  and  could  serve  as  the  basis  for 
a  more  complete  analysis. 

In  Sec.  IV  we  discussed  a  class  of  laser  guided  stationary  solutions  of 
both  the  integral  equation  formulation  and  the  two-dimensional  formulation  of 
a-c  boundary  motion.  The  problem  remains  of  Identifying  spatial  and  temporal 
Initial  conditions  which  allow  the  boundary  motion  to  settle  Into  such  a 
constant  velocity  stationary  state.  The  solution  of  this  problem  could  have 
Important  Implications  for  the  possibility  of  producing  uniformly  aligned, 
laser  crystallized  semiconductor  films.  In  studying  solutions  of  the  Integral 
Eq.  (3),  as  discussed  In  Sec.  IIA,  we  noted  that  for  a  given  Interval  fit  used 
In  numerical  Integration,  values  of  normalized  laser  scan  velocity  V  could  be 
found  which  were  large  enough  to  generate  damped  oscillatory  behavior  In  S(t). 
After  a  few  cycles,  S(t)  settled  down  to  a  constant  value,  consistent  with  the 
a-c  boundary  riding  at  a  fixed  distance  In  front  of  the  scanning  laser  Image, 
as  expected  In  a  laser-guided  stationary  solution.  However,  If  V  was  now  held 
fixed  and  5t  decreased  further,  the  damped  oscillatory  behavior  was  again 
replaced  by  perodlc  relaxation  oscillation  behavior  of  the  type  shown  In  Fig. 

2  for  n  =  0.3  or  0.6.  These  numerical  results  are  to  be  expected,  since  a 
given  value  of  6t  gives  rise  to  an  effective  finite  Initial  velocity  Veff  of 
the  a-c  boundary,  and  for  a  value  of  laser  velocity  V  >  Veff,  damped 
oscillatory  motion  should  result.  On  the  other  hand,  Eq.  (3)  Is 
characterized  by  singular  behavior  of  the  Initial  velocity,  so  that  decreasing 
5t  sufficiently  will  Increase  Veff  to  the  point  where  Veff  >  V,  and  the  a-c 
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boundary  again  outruns  the  laser,  resulting  In  periodic  relaxation 
oscillations.  As  we  have  seen,  this  singular  behavior  of  the  Initial  velocity 
of  the  a-c  boundary  Is  unphysical,  and  would  probably  be  eliminated  from  the 
mathematical  description  by  replacing  n  In  Eq.  (3)  by  a  velocity  dependent 
rigff  given  by  Eq.  (22).  Under  these  circumstances,  the  nature  of  the  solution 
obtained  should  not  depend  on  the  Interval  St  chosen  for  numerical 
Integration,  provided  that  St  was  sufficiently  small.  The  argument  Implies 
that  for  a  small  value  of  laser  scan  velocity  V,  the  Initial  motion  of  the  a-c 
boundary  would  be  characterized  by  some  finite  velocity  V0.  For  laser  scan 
velocity  V  >  V0,  we  would  then  expect  a  damped  oscillatory  behavior  of  S(t), 
with  the  motion  settling  Into  a  laser-guided  stationary  state.  A  slow 
decrease  In  V  could  perhaps  then  lead  the  a-c  boundary  motion  Into  a 
laser-guided  stationary  state  at  a  lower  velocity,  with  the  a-c  boundary 
riding  close  to  the  thermal  Image  of  the  scanning  laser.  However,  the 
stability  of  such  states  remains  questionable.  All  of  these  points  are  highly 
speculative,  and  require  further  exploration. 
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APPENDIX  A  -  Numerical  Solution  of  Integral  Equation 

The  non-linear  Integral  equation  (3)  cannot  In  general  be  solved 
exactly,  and  approximate  solutions  must  be  obtained  numerically.  Although 
the  results  are  not  entirely  satisfactory,  to  solve  Eq.  (3)  we  have  used  a 
slightly  modified  rectangular  rule  to  evaluate  the  Integral,  making  use  of  the 
already  determined  values  of  S( t")  at  earlier  times  x'  to  evaluate  S(t  ) 
self-conslstently  at  time  t.  The  modification  of  the  rule  consists  of  removing 
the  singularity  In  the  Integrand  at  x"«  x  by  treating  all  terms  except 

in  the  integrand  as  constant  In  an  Interval  fix',  and  Integrating 
(x-x')l/2  over  the  Interval  6x'.  The  Integral  equation  (3)  is  then  replaced 
by  the  approximate  equation. 


1 


a  exp 


n-1 

l 


m=0 


fnm 


S  +  V 
m 


exp 


(Tn-tm)1/2 


(Al) 


e*P  ]"[Sn_sm  +  v  (Tn_1m)F  /(‘Ih-''m)  !  » 


where  ynm  =  2(n-m)1/2  [(n-m)l/2  -  (n_m-i)l/2  ]f  and 

=  ^m  -  Sm-l  ,  m  *  0;  =0. 

6x 

Sm  Is  the  value  of  S  at  t  «  %,  and  (Al)  Is  solved  for  Sn  using  the  values  of 
Sn,  determined  for  earlier  times.  When  §m  <  -V,  the  assumption  Is  made  that 
the  front  has  ceased  to  move,  and  (§„,  +  V)  Is  set  equal  to  zero. 
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1 


Once  the  set  of  values  So,  Si,  ...  has  been  determined,  the  reduced 
temperature  at  point  S'  at  time  t  ■  tn,  given  by  Eq.  (5),  can  be  determined 
numerically  using  an  approximate  trapezoidal  rule  which  yields  the  relation. 


2 

n-1 

§  +  V 

T'  (S',  t  )  -  a  exp  {-S'  }  +  n 

l  y  _ 

.  m 

n  m*0  nm  1/2 

(t  -  x  ) 


n  m 

exp  {-r(t  -  t  )}  x 
n  m 


exp  -  [  S'-  S*  +  V(n,  -  ■*)]  /(*  -  ■*)  6t  .  (A2) 


Equation  (A2)  has  the  advantage  that  It  gives  reasonable  results  for  S'  well 
away  from  Sn,  and  agrees  exactly  with  (Al)  for  S'  *  Sn.  However,  It  gives 
rise  to  a  spurious  spike  near  S'  ■  S„.  The  origin  of  this  spike  can  be  seen 
by  comparing  the  second  exponential  In  the  Integrand  of  Eq.  (5)  with  the 
corresponding  exponential  In  the  summation  of  (A2).  For  S'  vS  In  Eq.  (5), 
the  exponential  term  ♦  -•  as  t'  ♦  t  yielding  a  vanishing  contribution  to  the 
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Integral.  The  term  corresponding  to  t'  ♦  t  In  (A2)  Is 
2 

-[S'  -S  +  V  6t] 

_ w-1  _  .  yielding  a  non-vanishing  contribution  of  the  summation 

which  peaks  up  when  S'  ~S  -  V6t  and  gives  rise  to  a  spike  of  width 


AS'  ~  (6t)  .  The  spikes  have  been  substracted  out  In  the  plots  of 

T'(S',  t)  vs  S'  shown  In  Figs  7,  8,  and  9. 


Appendix  B  -  Two-Dimensional  Stationary  Solution  of  a-c  Boundary  Problem 
We  consider  a  two-dimensional  heat  flow  problem  In  which  a  moving 
source  of  heat  scans  at  a  steady  velocity  v  along  the  surface  of  a  film  (1)  of 
thickness  bi  deposited  on  a  substrate  (2)  of  thickness  b2  (see  Fig.  Bl).  If 
the  properties  of  the  film  are  the  same  whether  It  Is  In  the  crystalline  or 
the  amorphous  state,  we  need  not  at  the  moment  specify  the  state  of  the  film. 
We  assume  that  the  moving  source  deposits  heat  only  In  the  film  and  uniformly 
In  the  x-dlrectlon.  The  heat  flow  equations  In  regions  (1)  and  (2)  are  then. 


where  G(y-vt)  Is  the  heat  energy  delivered  to  the  film  per  cm3  per  second. 

If  we  consider  only  stationary  solutions  In  a  frame  of  reference  moving 
with  velocity  v,  then  Tj  and  T2  are  dependent  on  time  only  by  virtue  of  their 
dependence  on  (y  -  vt).  In  this  moving  frame  of  reference  these  equations 
become, 

i!ll  +  £l  +  -0  (Bl) 

3x2  3u^  3u  *1 
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Fig.  B-l.  Diagram  for  the  calculation  of  the  exact  stationary  solution 
of  moving  a-c  boundary  problem. 
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32T2  +  32T2  +  C2P2V  3T2  ,  0  (B2) 

3x2  au2  *2  3,1 

where  u  *  y  -  vt. 

The  heat  flow  out  of  the  film  surface  Is  zero: 

«  0  (B3) 

x=0 

The  boundary  conditions  at  the  Interface  between  film  and  substrate  require 
that  the  temperature  and  heat  flow  across  the  boundary  are  continuous: 


At  the  base  of  the  substrate  the  temperature  Is  T^: 

[t2[x,u)]  -  Tb  .  (B6) 

x*(bi  +  b2) 

We  now  assume  the  film  Is  so  thin  that  the  temperature  variation  across 
It  Is  quite  small.  Integrating  (Bl)  with  respect  to  x,  and  using  (B3)  and 
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(B5)  to  eliminate  aTi/3x  yields, 

K2  [aT2  1  +  +  PlCjV  aTi  +  G(u)  -  0  (B7) 

biKi  "1x1  3u2  Ki  3u  Ki 

x*bj 

Finally,  If  we  allow  bi  to  become  very  small  compared  to  b2,  we  can  write  the 
coupled  equations  to  be  solved  for  Ti  (u)  and  T2  (x,u)  as. 


3^t2  +  fit  +  v_  £l  .  0  (B8) 

3x2  au2  *2  au 

*2  +  JL  li  +  Hal  -  0  (B9) 

bjKi  I  3x  Jx*0  au2  iq  au  *1 

with  boundary  conditions 

Ti(u)  =■  [T2(x,u)]  (BIO) 

x»0 

CT2(x,u)]  »  Tb  .  (Bll) 

x*b2 

We  Introduce  the  Fourier  Integrals  of  (Ti(u)-Tb)  and  (T2(x,u)  -  Tb): 
ti(k)  -  I"  (Ti(u)  -  Tb)  e<ku  du  (B12) 

-  m 

T2(x.k)  -  /"  (T2(x,u)  -  Tb)e^u  du  ,  (B13) 


58 


with  Inverses 


Tl(u)-Tb  -  1  /  nOO  e-*ku  dk  , 

77"  -“ 


T2(x,u)-Tb  »  1  /  t2(x,k)  e-*ku  dk 

TiT  -“ 


Applying  (B12)  -  (B15)  to  (B8)  -  (Bll),  we  obtain 


(B14) 

(B15) 


3*2  <^j  t2  =  o 

3x2 


(B16) 


*2 

blKi 


3*2 


3X 


‘“I  '1  *  9 

0  KJ 


T1  =  Ct2] 


[t2]  .  0 

*>2 

where  g(k)  is  the  Fourier  transform  of  G(u), 
g(k)  *  /  G(u)  e1ku  du  , 

.  oo 

a?  *  k2  +  1  vk  ,  and 
2  * 

oj2  >  k2  +  1  vk  . 

*1 


(B17) 

(B18) 

(B19) 


(B20) 
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The  solution  of  (B16)  Is  of  the  form, 

t2  =  Ae<«x  +  Be- °CX  .  (B21) 

Substitution  of  (B21)  Into  (B18)  and  (B19)  permits  the  expression  of  A 
and  B  In  terms  of  q.  Further  substitution  of  these  expressions  Into  (B17) 
gives  the  solution  for  q  In  the  form, 

T  * _ 9/Kl _  (B22) 

o£  +  a  K2  coth  (a,b-) 

1  2  biKi  ^  2 


The  Inverse  Fourier  transform  (B14)  leads  to  the  result, 

_ .  >  T  ,  r  gw/Ki  e"1ku  dk 

Tl(u)-Tb+1  /  1  _  ,  , 

dTf  -•  of  +  a  *2 coth  (ogb2) 

1  2  Wi 

He  will  take  G(y-vt)  *  G(u)  to  be  of  the  form, 

G(y-vt)  <*  G  •  9 (u)  where  0(u)  *  0,  u  <  -D,  u  >  D, 

e(u)  *  1,  -D  y  ••  <  D  . 

This  form  yields, 

g(k)  *  2  sin  kD  «  G 
k 


(B23) 


(B24) 


(B25) 
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It  Is  Important  to  note  that  for  several  heat  sources,  g(k)  Is  the 
superposition  of  contributions  due  to  the  Individual  sources.  Thus,  for  a 
case  In  which  separate  laser  and  latent  heat  Inputs  are  present,  the  solution 
of  the  problem  Is  the  sum  of  Individual  terms  of  the  form  (B23). 

With  g(k)  taken  to  be  of  the  form  (B25),  we  will  write  the  Integral 
(B23)  In  terms  of  normalized  quantities.  We  Introduce  a  characteristic 
normalization  length  d  to  be  specified  and  write, 

00 

Ti(S-)  =  Tb  +  lp  /  (Gd2/Ki)  2  sin  qD'  e-^  dq  (B26) 

TiT--  q  -4r(q)  +  0Q(q)  cot  Q(q) 

where  q  *  kd,  D'  -  0/d,  S'  -  u/d, 

21r(q)1/2  ■  oqd  ■  [q2  +  41Vq]J/2, 

Q(q)  =  -1(<*d)  (b2/d)  -  2  [r  +  2V  {M-l }  {-V  t  [p  +  V2  ]l/2 }  ]l/2  (b2/d),  (B26a) 

and 

B  ■  W2  . 

Klblb2 

M  =  *l/«2  , 

V  =  dv/4tq  . 


The  quantities  B,M,V  are  the  Important  parameters  of  the  model.  For 
convenience  we  will  In  this  appendix,  choose  d  *  b2  so  that  &  becomes 


K2b2/Klb1*  Tbe  (*)  s^9n  1n  Q(<l)  represents  the  fact  that  for  a  given  value  of 
r(q),  there  are  two  values  of  q  given  by, 
q  =  21  {-V  ±  [r  +  V2  ]l/2 > 

with  the  +  sign  corresponding  to  postlve  Imaginary  q  and  the  minus  sign 
corresponding  to  negative  Imaginary  q. 

m 

The  notation  p  /  signifies  the  principal  part  of  the  Integral: 

00  -  £  00 

p  /  dq  =  /  dq  +  /  dq,  where  e  Is  an  Infinitesimal  quantity.  This  Is 

-OO  -OB  +e 

Important  In  evaluating  (B26),  since  the  Integral  Is  singular  at  q  *  0. 

The  Integrand  In  (B26)  Is  a  single-valued  function  of  q  with  poles 
given  by  q  -  0  and 

D(q)  =  -4  r(q)  +  0Q(q)  cot  Q(q)  =■  0  .  (B26b) 

If  Q(q)  Is  complex,  as  It  will  be  for  much  of  the  complex  q-plane,  Q(q)  cot  Q(q) 
Is  still  well  defined.  For  Q(q)  pure  Imaginary,  QcotQ  ♦  |Q|  coth  |Q|. 

The  Integral  can  be  evaluated  by  contour  Integration  after 
appropriately  closing  contours  at  |q|  *  «.  Sin  qD'  Is  written  as 
(e1qD'  _  e-1qD ")/21 ,  and  the  Individual  Integrals  are  of  the  form 


p  /  e-1(lz  f(q)  dq  . 

For  Z  >  0,  the  contour  Is  closed  In  the  lower  half  plane,  while  for  Z  <  0,  the 
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contour  is  closed  In  the  upper  half  plane.  The  Integrals  can  be  evaluated 
using  the  residue  theorem.  The  contributions  from  the  zeros  of  D(q)  are 
obtained  by  writing  D(q)  In  the  vicinity  of  the  n-th  pole  as 

D(q)  2  ®iil]  (q  -  q  )  , 

n  aq  q  =  q  n 

L  J  n 

which  Introduces  a  factor  of  ( 30(q)/3q)q*q  Into  the  denominator  of  the 
expression  for  the  n-th  residue.  For  large  values  of  the  parameters  B,  the 
case  of  Interest  here,  the  zeros  of  D{q)  lie  on  the  Imaginary  q  axis  in  both 
the  upper  and  lower  half  plane  (corresponding  to  r(q)  real). 

Without  going  through  all  of  the  arithmetic,  we  present  the  results  of  the 
evaluation  of  (B26)  after  performing  the  contour  Integrations.  Ti(S')  can  be 
written  as, 

T i (S ")  -  Tb  +  l  (Gd2/K1)  (Fn+  (S'  -  O')  -  Fn+(S'+D')}  (S'  <  -  O') 

n+ 

Tl(S')  «  Tb  +  l  (GdZ/K!)Fn+  (S'  -  O')  (B27) 

n+ 

-  I  (Gd2/K!)Fn_  (S'  +  D')  +  (Gd2/Ki)  (_D'  <  S'  <  D') 

n-  6 

Tl($')  -  Tb  +  l  (GdZ/Ki)  (Fn.  (S'  +  D')  -Fn  -(S'  -  D ') }  (D'  <  S') 

n- 
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<  < 


where 


e[-V  ±  (r„*  +  v2)l/2  ]2Z 

Fn*  (Z)  3 _ 

(30±/3r)n±  Or/3q)n±  21  [-V  t  (r  +  v2)  1/2] 

n* 


(B28) 


0  =  -4  r(q)  +  6  Q  (q)  cot  Q  (q)  , 

±  ±  ± 

-4r(q)  *  q2  +  4  IVq,  and 

Q  (q)  *  2[r  +  2V  [M  -  1]  {-¥  ±  ( r  +  v2)l/2 }]  1/2. 

♦ 

The  rn±  correspond  to  the  zeros  of  D±  .  The  values  of  q  corresponding  to  rn± 
are 

q„±  3  2  1  {-V  t  [rn±  ♦  v2]l/2}  . 

The  r's  are  real  numbers  while  the  q's  lie  along  the  Imaginary  axis. 

Me  will  be  concerned  with  the  case  where  D'  Is  very  small.  Expanding 
Fn±  (s'  -O')  and  Fn*  (S'  +  O')  to  first  order  In  D*,  and  taking  the  remaining 
derivatives  In  (B28),  we  obtain 


Tl(S') 


Th  +  l 

D'  (GdZ/Ki) 

X 

n+ 

[2  -  e 

cot  Q  2  1  ( 

V  -  CSC  Qn+ 

1  +  v  ! 

] 

PST"  j  | 

[  /(  rn+  +  v2)| 
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» 


e-2VS-  e-2  [rn+  ♦  V2]1/2  |S'| 

Crn+  +v2]1/2 


(S'  <  o)  , 


Ti(S') 


Tb  +  l 
n* 


0'  (GdZ/Ki) 

[2  -  6 

cot  Qn — csc2  Qn_ 

ll-  V  (M-l)  j 

On- 

1  yrrrTvJTj 

t 


e-2VS'  e-2  [ft-  ♦  V2  ]1/2  |S'| 

[rn.  +  v2]i/2 


(B29) 


(S'  >  0). 


Now,  suppose  that  the  source  of  heat  energy  delivered  to  the  film  Is 
the  latent  heat  given  up  as  the  a-c  boundary  moves.  The  energy  G  delivered 
per  cm^  per  second  Is, 


G 


Lpiv 

“Iff 


where  2D  Is  the  width  of  the  region  over  which  this  energy  Is  delivered  In 
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time  2D/v.  Introducing  normalized  quantities, 

(6d2/K1)  -  2  L  V  .  (B30) 

Ci  TP 

Substitution  of  (B30)  Into  (B29)  and  some  simple  manipulation  of  the 
resulting  expression  leads  to. 


Ti(S') 


Th  +  LV  v  e-2VS'  e-2[rn+  *  V2]1/2  |S'|  x 

CT  - - 

n+  [rn+  +  v2]l/2 


(B31a) 


_ 1 _ 

ti-*  js*$  .  ««*  v| 

(S-  <  0)  , 


MS')  «Th  +  LV  l  e-2VS *  e-2[rn-  +  V2]1/2  |S'|  (631b) 

1  b  ^  n-  [r~  +  V?ji72 
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**- 


1 


[1  -  e/2  | cot  Qp-  -  esc2  Qn“ 
Qn- 


,1  .  vfej) .  ] 


(s-  >  o). 


Equations  (B31a)and  (B31b)  are  used  to  calculate  the  velocity  of  the  a-c 
boundary  In  the  steady-state  runaway  condition  In  Sec  1 1  IB. 
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